Bosonic Reduction of Susy 
Generalized Harry Dym Equation 

Ashok Das° and Ziemowit Popowicz 6 

a Department of Physics and Astronomy 
University of Rochester 
Rochester, NY 14627 - 0171, USA 

b Institute of Theoretical Physics 
University of Wroclaw 
pi. M. Borna 9, 50 -205 Wroclaw, Poland 

February 8, 2008 
Abstract 

In this paper we construct the two component supersymmetric generalized 
Harry Dym equation which is integrable and study various properties of this 
model in the bosonic limit. In particular, in the bosonic limit we obtain a 
new integrable system which, under a hodograph transformation, reduces to 
a coupled three component system. We show how the Hamiltonian structure 
transforms under a hodograph transformation and study the properties of 
the model under a further reduction to a two component system. We find a 
third Hamiltonian structure for this system (which has been shown earlier to 
be a bi-Hamiltonian system) making this a genuinely tri-Hamiltonian system. 
The connection of this system to the modified dispersive water wave equation 
is clarified. We also study various properties in the dispersionless limit of our 
model. 
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1 Introduction 



The Harry Dym equation |2| is an important dynamical system which 
is integrable and finds applications in several physical systems. Together 
with the Hunter- Zheng equation [HI IH El El E|, they define a hierarchy of 
integrable systems for both positive and negative integers. This system has 
been vigorously studied in the past as well as more recently from various 
points of view. In particular, the supersymmetrization of such a system has 
been discused systematically in [311] and a two component generalized Harry 
Dym equation has been constructed in [TU] . 

The construction of supersymmetric integrable systems and the under- 
standing of their properties is important for a variety of reasons [TT] IT!?! E] ■ 
One of the interesting features lies in the fact that supersymmetric integrable 
models can lead to new integrable systems in the bosonic limit ^3JE1E3- As 
we have already pointed out in connection with the N = 2 supersymmetriza- 
tion of the Harry Dym equation, one of the supersymmetrizations leads to 
an interesting new coupled integrable system in the bosonic limit jHj- It is 
with this aim that we have chosen to construct and study the two component 
supersymmetric (susy) generalized Harry Dym equation in this paper. The 
construction is quite interesting and leads to many new features including 
the fact that in the bosonic limit we obtain new integrable equations. 

Our paper is organized as follows. In section 2, we recapitulate briefly 
the two component generalized Harry Dym equation and some of its prop- 
erties. In section 3, we construct the two component supersymmetric Harry 
Dym equation which is integrable. The Lax representation as well as the 
Hamiltonian structures are discussed in detail in terms of different variables. 
We construct the hodograph transformation for the bosonic limit of such a 
system. In section 4, we address the question of how a Hamiltonian struc- 
ture transforms under a hodograph transformation. Using this, we derive 
the Hamiltonian structure for the transformed equation and show that un- 
der some approximation, the system reduces to a three component coupled 
MKdV system of equations. A further reduction, in section 5 takes this equa- 
tion to one which has been studied earlier and we find a new Hamiltonian 
structure for this system which genuinely makes it a tri-Hamiltonian system. 
The connection of this system with the modified dispersive water wave equa- 
tion is also clarified. In section 6, the dispersionless limit of our system of 
equations is studied systematically and various associated properties includ- 
ing polynomial and nonpolynomial charges are derived. We conclude with a 
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brief summary in section 7. 



2 Harry Dym and the Generalized Harry Dym 
Equations: 

As is well known, the Harry Dym equation [1] can be written in the form 

dw 3 

— = w w xxx , (1) 

where subscripts denote derivatives with respect to the corresponding vari- 
ables and the equation can be obtained from the following Lax representation 

^ = 4[(L)>| ,L\, (2) 

with n — 1 and 

L = w 2 3 2 . (3) 

Here > 2 refers to the projection onto the sub-algebra of the pseudo differ- 
ential operator P, namely, 

(P)> a = 5>*. (4) 

i=2 

The hierarchy of Harry Dym equations is integrable and finds applications 
in various physical examples. 

The Harry Dym equation can be generalized to two dynamical variables 
in the following way [TO] ■ Let us consider a Lax operator in the product form 

L = w 2 d 2 u 2 d\ (5) 

where w, u denote two dynamical variables. It is straightforward to check 
that the non-standard Lax representation 

leads to the system of coupled equations for w and u of the form 

dw 
~dt 

f = « 3 (-- 1/2 » 3/2 )„- p> 

KJ v \ / XXX 



wHw- 1 ' 2 ^' 2 



This is the two component generalized Harry Dym equation which can be 
written in the Hamiltonian form as 



w \ i w 3 d 3 u 3 \ I §g 

3Q3 W 3 o I 1 SH 



(8) 



Su 



U I \ U O W 

where 

H = -2 J dx (wuj 5 . (9) 

We note that we can rewrite the system of equations (J7j) in a simpler form 
in terms of the new variables 

u = ae b , w = ae~ b , (10) 

as 

a t = a 3 (a xxx + I2a x b x + 126^^0), 

b t = -2a 2 ((ab) xxx - a xxx b + 4b 3 x a). (11) 

This equation can be writen in the Hamiltonian form 

a \ I / ar (X - X^) a 3 -a 3 (X + X^) cr 



,3 (y V\\ ^3 ^3 (v i Vt^ „2 \ f SH 

,2 ) 1 6H_ 

Sb 



(12) 



b J 4 ^ a 2 (X + Xt) a 3 -a 2 (X - X+) a 
where the operator X can be identified with 

X = e 2b d 3 e- 2b , (13) 
with X^ representing the Hermitian conjugate and 

H = -2 I dxa- 1 . (14) 



3 Supersymmetric Generalized Harry Dym 
Equation: 

The supersymmetric extensions of the Harry Dym equation have already 
been constructed in jH] . Let us briefly recapitulate some of the features of the 
N = 2 extended supersymmetric Harry Dym equation which is most useful 
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from the point of view of constructing the generalized supersymmetric Harry 
Dym equation. In this case the bosonic dynamical variable is generalized to 
a bosonic N = 2 superfield 

W = w + O 1 xi + 6 2 X2 + O 1 e 2 w u (15) 

where W\ is the original bosonic dynamical variable while wq is the new 
bosonic dynamical variable and Xi t i = 1 > 2 are the new fermionic dynamical 
variables necessary for N = 2 supersymmetry. We note that 9i,i = 1,2 
represent the two Grassmann coordinates of the N = 2 superspace. We can 
define the two supercovariant derivatives on this space as 

d d d d 

which satisfy 

{D u D 2 } = 0, D\ = D\ = d. (17) 

In |H] we showed that it is possible to construct four different N = 2 
supersymmetric Lax operators that lead to consistent equations. Namely, 
the general Lax operator with the N = 2 superfield T 

L = T- l d 2 + k l (D 1 T- l )d + k 1 (D 2 T- 1 )d + 

(k 2 (D 1 D 2 T)T~ 2 + k 3 (D 1 T)(D 2 T)T- 3 )D 1 D 2 , (18) 

leads to a consistent non-standard Lax representation 

dL r a i 

— =4[(L)| 2 ,L], (19) 

only for k\ = k 2 = k 3 = or k\ = k 2 = — -y = 1 or ki — |, k 2 — — 
or k\ = k 2 = \ik% = |. We will consider here only the last case which 
allows us to carry out the construction of the two component supersymmetric 
generalized Harry Dym equation. 

Let us note (which has been pointed out in [S] as well) that in this par- 
ticular case, the Lax operator (|18|) can be written as 

L = -(WD 1 D 2 ) 2 , (20) 

where W is related to T and the Lax equation (J19)) leads to the dynamical 
equations 

W t = i(6(D 1 D 2 H/)( J D 1 D 2 H/ x )H/ 2 - 2W XXX W 3 

+Z(D 1 W XX )(D 1 W)W 2 + 3(D 2 W XX )(D 2 W)W 2 ). (21) 
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This equation can be written in the Hamiltonian form 



*=*m- (22) 



where 



T>i = -^W 2 D 1 D 2 dW 2 , H = ~ J dxd 2 6 {D 2 W){D 1 W)W~\ (23) 

with d 2 6 = d9 1 d9 2 . 

In order to construct the two component supersymmetric generalized 
Harry Dym equation, let us consider the two differential operators on the 
N = 2 superspace 

L x = FD 1 D 2 , L 2 = GD 1 D 2 , (24) 

each of which has the form of the square root of (|2(J|) . Here F, G denote two 
bosonic superfields on the N = 2 superspace. Let us next construct a new 
Lax operator as the product 

L = -L X L 2 . (25) 

It can be checked that the two component susy generalized Harry Dym equa- 
tion follows from the non-standard Lax representation 

§ = 4[L| 2 ,L] (26) 

In fact, the equations are derived in a much simpler manner with a parame- 
terization of the form (jl(Jj) . namely, 

F = e w U, G = e~ w U. (27) 

With this parametrization, the dynamical equations can be written in the 
Hamiltonian form 



U \ 1 / U 2 (X - X f ) U 2 -U 2 (X + U 



8H 2 
5U 

W J 4 ^ U (X + Xt) U 2 -U (X - Xt) U J \ §k 
where X = e w DiD 2 de~ w and 

[DzU^DtUy 



(28) 



H 2 = TrfLM = / dxd 2 6 



(D 2 W)(D 1 W)U - 



U 



(29) 



It is also possible to write these equations in a still simpler form by intro- 
ducing the variables 

t/ = ^U W=Uhig-hif), (30) 

vfg 2 

in terms of which we have 

/ 

, 9 

where 

Y = D 1 D 2 d, H 2 = -4 [dxd 2 6-^= D X D 2 ( ^— 




(31) 



It is now obvious that the Hamiltonian structure in ()31|) satisfies Jacobi 
identity. We will, however, continue to work with the first representation of 
the Hamiltonian equations (|28|) . This system of equations is integrable and 
the conserved charges can be obtained from the traces of the Lax operator 
in the standard manner. 

In the bosonic limit where we can restrict the superfields and the Hamil- 
tonian to the forms 

W = w + 9 1 9 2 wi, U = u + 9 1 9 2 u 1 , 

H 2 = Jdx (wlu + wjuo - , (32) 

equation (|28j) leads to a system of four interacting equations of the form 

w ot = - (2woxxxul + 6w 0xx u 0x ul - w^ x ul + 3w 0x ul x u 
+3w Ox ulu - 3wlwo x ul) , 

wu = 2 {-6 w 0x X v>ixul - ^ w ox u i u o - 6w 0x u lxx ul - 12w Qx u lx u 0x uo 
+3w 0x ul - QwoxUiUoxxUq + 3w 0x uiul x + 2w lxxx ul 



+12w lxx UoxUo + ^w lx u 0xx u + 12wixU 0x u + Qw lx w 0x u 

2 

dxxx^Q 



-6w lx w 2 ul + Qwfuoxul - 3w 2 woxUiul + Qwiu 0x xxul 



+ 12w!UoxxUOxUo - QWiUix^Uo + QWiWQxxWoxul + ^WiW^UoxUl) 

uot = {uoxxxul - 3u lx u x ul + 3w lx w 1 u [ + 3wlu Qx ul) , 
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uu = - (2u lxxx ul + 6u lxx u 0x ul + 6u lx u 0xx ul - l2u lx u\u Q 

+Quiu 0xxx ul + Qwo^wo^ul + Qwl x u lx ul + Qw^UxUqxuI 
-6w 1xx w 0x Uq - Qwi x w 0xx Uq - 24w lx w 0x u 0x ul + \2wi x wiUiul 
-3wfw 0x v,Q + I2wluxu 0x ul - §WiW Qxx u 0x ul - 3wiwI x Uq 
-12wxWq x uq xx uI - §wiw 0x ul x ul + 3wxw 0x ulul) . (33) 

When Wq — W\ — 0, this system of equations reduces to the bosonic limit of 
the N = 2 susy Harry Dym equation (|2Tj) . On the other hand, when we set 
w = Ui = 0, and identify u = a, Wi = 2b x the system of equations goes over 
to the generalized Harry Dym equation (jllj) (except for an overall factor in 
the second equation which we are unable to get rid of by any scaling) . 

Let us next study the behavior of this system of equations under a hodo- 
graph transformation jT^] of the variables (x,t) to (y, r) defined through 



x = Po(y,r), r = t, (34) 



so that we have 



d dy d Id 

dx dx dy po y dy ' 

— - — + — = — + ^L— (35) 

dr dt T dx dt p 0y dy ' 

Defining the variables 

= 9o ; w i = Qi, u o = Po y , Ui = pi, (36) 
the system of equations (|3^|) goes over to 

POr = g i 2 P0yyy ~ ^lyyPoy ~ MPOy + ZQlPly) , 

Pit = 2 { 2 Piyyy + 6 PiyPoyyyPoy ~ ^PiyPlyyPvy ~ ^PiyPl 

+QPlP0yyyyP y ~ 24:p l p 0yyy P0yyP y 2 + ^PlPlyyPa'y + QqOyyQOyPl 
+6<l0yPly ~~ ^QlyyQOyPOy ~ ^QlyQOyyPOy ~ ^-2qi y qQ y p 0yy 



+ I2q ly qip 1 p 2 0y - 3qlq 0y ply + 3qlp ly pl y + 12gj>ip 0m Poy 

SqiqOyyPOyy ~ Ml^yPOy ~ ^QlQOyPOyyy + ^QlQOyPoyyPoy 
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+3q 1 q 0y p 2 1 p y) , 



<?0t 



9l7 



+3q yPlPoy - SqOyPlPOyyyPoy + $%yPlPlyyPoy + ^Qlyyy 
+6qiyyP0yyP0y + 12 QlyPOyyyPoy ~ ^QlyPlyyPoy ~ ^QlyPl 
+6<lly<l0y + ^QlyQlPly + GqfPOyyPOy ~ ^QOyPlPOy 
+6<llPOy y yyPoy - \1qxPVyyyPVyyPvy + ^iPlyyPoy ~ ^QlPlyPl 

+6qiqo yy qo y + Qqiqo y Po y yPoy) ■ 



(37) 



We note from (|3*7)l that the equation for the variable q is decoupled. 
As a result, we can set go — for simplicity. The Hamiltonian structure 
for the system of equation involving wi,Uq,ui can be obtained from the 
Dirac reduction of the bosonic limit of the Hamiltonian operator in (|28|). In 
the simple case of two variables, for example, the Dirac reduction works as 
follows. Let us assume the equations of motion to have the form 



u 



P P 

1 uu 1 u 

p p 

mi. iv 




(38) 



a a j Puvi Pvui Pw repre 



where u, v denote the two dynamical variables with P u 
senting the elements of the Hamiltonian operator. If we can set v = 0, then 
it can be verified directly that j2j E| the Hamiltonian structure reduces to 



u t 



P - P P _1 P 

1 UU UV 1)1) 1 vu 



5H 

5u 



v=0 



for the reduced system. 

For the three component system under study in (J3~3~|) (with w 
Dirac reduction leads to 







( 


SH 2 
Suo 


\ 


Ui 


= V 




5H 2 
Sui 




\w 1 J 


t 


\ 


&H 2 
5wi 


J 



(39) 



0), the 



(40) 



where the matrix Hamiltonian operator T> has the elements 

pa. 1 ) = du 4 + 4d, 
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p(l,2) = 4 UlyU 3 _|_ 4 WlM%M 2 _|_ 4 MlM 3^ ) 

D^' 3 ) = 2u; ly no + 2w 1 u 0y ul + 2wiuld, 

p(2,2) = _ a 3 w 4 _ ^3 _ 0(12^^ _ 4m 2 m 2) _ _ 4^2)^ 

p(2,3) _ 4 u , 1 ^ MlM 2 _|_ 4 U ; lMlUo ^ Uo _|_ 4^iniWQC?, 

p( 3 > 3 ) = d 2 ul + u 2 d 3 + d(w 2 1 ul-3u 2 0y ) + (w 2 1 u 2 ~3u 2 0y )d, (41) 



and 



tf 2 = J dx lw 2 lUo - ^Ij . (42) 

4 Transformation of the Hamiltonian Struc- 
ture under the Hodograph Transformation: 

In this section, we will discuss how the Hamiltonian structure transforms 
under a hodograph transformation. Namely, we already have the Hamilto- 
nian structure for the three component equation in (JHHj) with w = 0. We 
would like to determine the Hamiltonian structure for the three component 
equation in ()37j) with q = which is obtained from under a hodograph 
transformation. We note that the Hamiltonian (J42j) in terms of the three 
variables po,Pi and q\ takes the form 

H = S dV {- V -pf y - pl + ^- (43) 

Furthermore, from the definition of the hodograph transformations in (J33j) 
and (jnnj), we obtain 

f u o \ ( Po \ 

where S has the form 



(45) 



P0y9~ 


y 





o\ 


Plyd~ 


POy 


1 





Qi y d~ 







1 ) 
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The gradient of the Hamiltonian transforms under the hodograph transfor- 
mation as 



(46) 



8H 

Su 






f 


8H 
Spo 




SH 




= K 




SH 




<5tii 






Spi 




ML 


j 




\ 


5H 


J 


Swi 




Sqi 



where 



K 



V 








POy 





It follows now that 

/ Po \ ' m 

Pi 

\ 11 ) T 



/ SH 2 \ 

8uq 



sv 



8H 2 
5u i 
SH 2 



SVK 



V Sw\ ' 





Poy 

/ SH 2 \ 

Spo 
5H 2 
<5pi 

\ % / 



Po, 



(47) 



/ SH 2 \ 

<5po 
SH 2 
Spi 

\ SJ ^ I 

V 591 / 



(48) 



The form of the transformed Hamiltonian structure X> can now be determined 
easily to have the following elements 



£,(1,1) 

p(1.2) 
p(l>3) 
p(2,2) 



-2pQ y d~ l q X y + 2gip 0y , 

- 2 ( d PoyyyPoy + Po yy yP y d ) + §(. d Pl yy Poy + PlyyPQy d ) 

-2 Ply d- l Ply + A{dp\ + p 2 <9) - 2«9 3 , 
-2p ly d~ x q lv + 2gipi y + 4ftj?i0), 

-2q ly d- l q ly + d 2 p^d + a Po " 2 a 2 + <9g 2 + qfd. (49) 
We can further simplify the equations (j37j) defining the new variables 

POyy 



p(3,3) 



/ 



9 = Pu 



QlPOy, 



(50) 



which leads to the equations 

/t = (2f yy -f-Qg y g-3g 2 f + Qh y h + 3h 2 f 
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9r = fayy ~ % 3 + Sgfy - 3fl-/ 2 + 3h 2 g J + 6/i y %, 

/i r = (2h yy + 5h 3 + 6hf y - 3hf 2 - 3hg 2 ) -6hg y g. (51) 
\ 'y 



The Hamiltonian structure D, in this case, transforms to 





-2/A 


- 1 f y + df + f 2 d-2d 3 , 




p(l,2) = 




1 g y + 2d(dg + gd) + 2dgf + 2gdf, 




P(L3) = 


-2/„a- 


- l h y + 2d(dh + hd) + 2dhf + 2hdf, 




p(2,2) _ 


-2g y d- 


l g y + d(f 2 + V - 2f y ) + (f 2 + V - 


2/,)9-29 3 , 


p(2,3) _ 


-2g y d' 


" 1 /i J/ + 4(9% + /i^), 


G 


p(3,3) = 


-2h y d' 


- 1 / iy + a(4/ i 2 -/ 2 + 2/ y ) + (4/ i 2 -/ 2 H 


-2/,)a + 2d 3 , 



and the Hamiltonian in the new variables has the form 

H = Jdy (-f 2 -g 2 + h 2 ). 

All the operators d in the expression (J52)) refer to derivative operators with 
respect to the variable y. We note that H = J dy f is a conserved quantity 
which defines the Casimir of the Hamiltonian structure T> in ()52j) in the sense 
that it annihilates the gradient of the Hamiltonian [S] . 

5 Connection with the Modified Dispersive 
Water Wave Equation: 

It is well known that the Harry Dym equation can be transformed to the 
MKdV equation under a hodograph transformation jTB]- Similarly, the gen- 
eralized Harry Dym equation can be transformed to two coupled MKdV sys- 
tems under a hodograph transformation ^H]- Such two component coupled 
MKdV systems have already been classified. In this section, we will study 
the connection of our system of equations with other equations. Let us note 
that if we set h — in (|37)l (which would correspond to setting W = for the 
super symmetric generalized Harry Dym equation) the system of equations 
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takes the form 

fr = (2f yy -f-6g y g-3g 2 f^ 

9r = (2g yy -5g 3 + 6gf y -3gf 2 ^. (53) 

This is a system of two coupled MKdV equations and under this reduction 
the Hamiltonian structure can be obtained from the Dirac reduction of the 
operator V as discussed earlier and we obtain the two by two matrix V with 
elements 

VQ-U = -2f y d- 1 f y + df + fd-2d 3 , 

= -2f y d- 1 g y + 2d(dg + gd) + 2dgf + 2gdf, (54) 
= - 2 g y d- 1 g y + d(f 2 + 4g 2 -2f y ) + (f + 4g 2 -2f y )d-2d 3 . 

We note that, as in the previous case, J dy f is conserved and defines the 
Casimir of the Hamiltonian structure. However, in the present case, H = 
J dy g is also conserved and can be used to construct a new system of equa- 
tions 




Under a change of variables 

f=(a + b), g = i(a-b), (56) 

the system of equations (|53*|l takes the form 

a T = (a yy + 3aa y — 3ba y + a 3 — 6a 2 b + 3ab 2 ) y , 

b T = (b yy + 3bb y — 3ab y + b 3 — 6b 2 a + 3ba 2 ) y , (57) 

which has been studied earlier by Sakovich JH] and Foursov JH]- On the 
other hand, under this change of variables (}56|) . the new system of equations 
(|55j) has the form 

a T = 2i(a y + a 2 — 2ab) y) 

b T = 2i(-b y - b 2 + 2ab) y . (58) 
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Theses constitute the system of coupled Burgers equations [2] and in the 
following we will neglect the overall factor 2i. 

Both of the system of equations ()57|) and ()58|) are Hamiltonian with the 
Hamiltonian structure (it is the structure (p>4*)) in the new variables) 

V { q ,1) = -2a y d~ l a y + d{a y + 2a 2 - 2ab) + {a y + 2a 2 - 2ab)d, 
pj 1,2) = -2a y d- x b y -d 3 -d(da - 2ab + a 2 + ad) + 

(db + 2ab -b 2 + bd)d + 2bda, 
Vf 2) = -2b y d- x b y + d(b y + 2b 2 - 2ab) + [b y + 2b 2 - 2ab)d. (59) 

It is already known that the systems of equations (joT|) and (Jo"H|) are bi- 
Hamiltonian [211^] with Hamiltonian structures T>i,T>2 which we describe 
below. However we find that, in fact, both these systems of equations are 
tri-Hamiltonian systems much like the two boson equation. For example, we 
can write the system of equations (Jo"%J) as 



a 



SHq \ / &ffi \ / 5H 2 

M £ \=Vi\ tl \= V A 1 I > (6°) 



Sb / \ Sb / \ 5b 



where 



H = J dy (a — b), Hi = J dyab, H 2 = J dy (a 2 b — ab 2 + a y b) . 

The other two Hamiltonian structures have been studied earlier to have the 
forms 



—2ad — a y d 2 + (a — b)d + a y 

-d 2 + (a - b)d - b y 2bd + b y 



, d , 

*=l 80 l- («) 

Similarly, we can write the system of equations (foTj) as 



a 



SHi \ _ / SH2 \ / 5H 3 

^ol 1 =VA * =V 2 £ 1, (62) 



6/ ffli \ 6H 3 

Sb / \ Sb / V 56 
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where 

6. = \ J W** - 3»\ - + 2 b a„ + Ua y a + 2ta-). (63) 

A careful analysis shows that the new Hamiltonian operator T> can, in fact, 
be written as 

V = -V{D 2 X V X . (64) 

This shows that the three Hamiltonian structures are related through the 
recursion operator 

K = Z>iX>2 \ 

as 

v x = nv 2 , v = -n 2 v 2 . 

Consequently, these define compatible Hamiltonian structures and the asso- 
ciated Nijenhuis torsion tensor vanishes [2D] . 

To close this section we point out that the system of equations (I58]l can 
be mapped to the modified dispersive water wave equations [2j under the 
transformation 

f 

a — > —v, b —> w — v, t — > — — . (65) 
Namely, under this transformation ()58|) goes into 

v T = -(-v y + 2vw-v 2 ) y , 

w T = ^{w x - 2v y - 2v 2 + 2vw + w 2 ) y , (66) 
which corresponds to the modified dispersive water wave equation. 



6 Dispersionless Limit: 

Given a dispersive system, one can obtain the dispersionless limit by taking 
the long wavelength limit [22] . In this limit the three component system 
obtained in (jST]) takes the form 

fr = (-f-Zg 2 f + Zh 2 f) y , 

g T = ^5g 3 -3gf 2 + 3h 2 g^ + 6h y hg, 

h T = ( + 5h 3 - 3hf 2 - 3hg 2 ^) -6hg y g. (67) 



15 



This can be written in the Hamiltonian form with the dispersionless limit of 
the Hamiltonian structure D in (|52j) 



£,(1,1) = 


-2fyd~ 


-'fy + df + fd, 


p(l,2) = 


-Vyd- 


1 g y + 2dgf + 2gdf, 


p(!' 3 ) = 


-2f y d- 


- x h y + 2dhf + 2hdf, 


X>( 2 > 2 ) = 


-2g y d' 


l g y + d(f 2 + Ag 2 ) + (f 2 + Ag 2 )d, 


p(2,3) _ 


-2g y d' 


^hy + Aidhg + hgd), 


p(3,3) _ 


-2h y d~ 


- 1 / ij/ + 9(4/i 2 -/ 2 ) + (4/ i 2 -/ 2 )a 



then in these variables the dynamical equations can be written as 



The Hamiltonian structure correspondingly takes the form 



p(l,l) = a ( 4F 2 _ _ 15QF j + ( 4F 2 _ _ 15GF }Q _ 2F y d~ 1 F y , 

= d{3G 2 + 20FG) + {3G 2 + 20FG)d + AFGy + 3G y G-2Fyd- 1 G 

p(i,3) = 208HF + (20HF + QHG)d + AFH y + 9HGy-2Fyd- 1 H y , 

V {2 ' 2) = d{AG 2 -16GF) + {AG 2 -lQGF)d-2G y d- 1 Gy, 

V {2 ' 3) = ~1Q(8HF + HFd) -8HGd-AHGy + 12H y Gy-2Gyd- 1 Hy, 

V {3 ' 3) = d(lQHF + Y2EG + 16# 2 ) + (1QHF + Y2EG + 16# 2 )<9 



and the Hamiltonian H = J dy (— f 2 — g 2 + h 2 ). 
If we now define the new variables 



F=h 2 - jG, G = -3g 2 + H, H — 3h 2 



(69) 




(70) 



— 2H v d H„. 



(71) 
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As we see, the first two equations do not depend on H and the third equation 
allows us to set H = consistently. In this case, with a redefinition of 
variables / = 3a — 2b and g = 4(6 — a), the first two equations reduce to a 
much simpler system, namely, 



5a 2 — Aab 



y 

b T = (-26 2 + 2ao) + 2ba y . (72) 

In this case, the Dirac reduction of the Hamiltonian structure becomes sin- 
gular. Nevertheless, we have explicitly verified, using the method of prolon- 
gation that the Hamiltonian structure 

h\P> _ Or, Pt^r, 

Vi 



T>^' 1] = 4<9(4a 2 - 3a6) + 4(4a 2 - 3ab)d - 2a y <T 1 a 



V { ^ 2) = Ad(2ab-b 2 )+A{2ab-b 2 )d-Aba y -Abyb-2a y d- l b 



Vq = 4(<90 2 + b 2 8) - 2byd- l by, (73) 

satisfies Jacobi identity and generates the dynamical equations as a Hamil- 
tonian system with the Hamiltonian H Q = | J r dy a. Furthermore, this new 
system of equations (J72|) is, in fact, bi-Hamiltonian 



V ° /a = V A /* > ( 74 ) 










5a 




8a 


SHo 




8b / 




\ 8b 



where H\ = J dy (20a 2 — 16a6) and 

~ J_ / da + ad db + bd \ . . 

1 ~ 12 V db + bd db + bd J ■ [ } 

We can now construct the recursion operator R = V~ 1 V and obtain the 
infinite series of conserved charges. The first few have the explicit forms 

H 2 = J dya(-2la 2 -8b 2 + 28a&), 
H 3 = J dy a(-429a 3 + 729a 2 £> - 432a& 2 + 64b 3 ) , 
Hi = J dy a(-2431a 4 + 5720a 3 6 - 4576aV + 1408ao 3 - 1286 4 ) . (76) 
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In addition to these polynomial conserved charges, we can also construct a 
nonpolynomial series of conserved charges using known techniques in hydro- 
dynamics, which we explain below. 
Let us note that 



#x = J dyJ^ - 1, H 2 = J dyaV^b (77) 

also define conserved quantities for the system. These charges have been 
constructed using the analog of the Tricomi equation used in the theory of 
polytropic gas dynamics [21] . Namely, if a conserved density H (H = J dy H) 
depends on a and b only, it satisfies H T = G y where G depends only on a 
and b. Using (J72*j) . this can be shown to lead to the following analog of the 
Tricomi equation 

AaH a>b + 2aH aa + 2bH bb + H b = 0, (78) 

where the subscripts denote derivatives with respect to the particular vari- 
able. The conserved quantities Hi in (J77)l can be obtained as particular 
solutions of the Tricomi equation. We note that if we scale H — ► (a — b)H 
then equation (J78j) transforms to 

4aif a>6 + 2aH aa + 2bH bb + 3H b = 0. (79) 

The general solutions of equations (|75|) and (fTTJj) can be written as 

71-1 

H n = Y,KnO n ~ k b\ (80) 

fc=0 

where 

X - 9\ (n — k + l)(n — k) 

" - 2Xk -^4k(n-k) + 2k(k-l) + zk (81) 

for k > with Ao, n = and z = 1 for (|78j) while z = 3 for ()79j1 . In the first 
case, we have polynomial charges the second leads to nonpolynomial charges 
(because of the factor \/a — b) . 

If we further change the variable b to b = 3c 2 , then equation (|72|) has the 
conservative form 

5a 2 - 12ac 2 



Ac A + 2ac) , (82) 
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which can be written in the bi-Hamiltonian form 



a 
c 



SH \ ( 5H 1 



Vo \ * )=V 1 [ s % ) , (83) 



5b 1 \ 5b 



where 

* = T* ( da t ad fa I ■ M 

while Po has the matrix elements 



V£' 1} = 4d(4a 2 -9ac 2 )+4(4a 2 -9ac 2 )d-2a y d~ 1 a. 



V^' 2) = -4(3c 3 -2ca)d + 2ca y -2a y d~ 1 c y , 
V (2,2) = dc 2 + c 2 d-2c y d~ 1 c y , (85) 



and H Q =lfdya,H 1 = j dy (20a 2 - 48ac 2 ). 

Finally let us comment on the different possible reductions of the system 
of equations (J57J) . We note that when h — 0, we have 

/r = (-/ 3 -3a 2 /)^, 

^ = (-5o 3 -3o/ 2 )^, (86) 

and it corresponds to the dispersionless limit of equation (|53|). We will not 
discuss this further. If we set / = 0, then leads to F — |(g 2 — /i 2 ), G = 
— 3(g 2 — /i 2 ), = 3h 2 and we have the system of equations 

G T = --(G 2 ) y , 



4 
1 
2 



H T = —-(HG) y + HG y , (87) 



where the first equation is decoupled. 



7 Conclusion: 



In this paper, we have constructed the two component supersymmetric gen- 
eralized Harry Dym equation which is integrable and have studied various 
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properties of this model in the bosonic limit. In particular, we find a new 
integrable model in the bosonic limit which under a hodograph transforma- 
tion maps to a system of three coupled MKdV equations. We have shown 
how the Hamiltonian structure transforms under a hodograph transformation 
and studied the properties of the system under a further reduction to a two 
component system. We find a third Hamiltonian structure for this system 
making this a genuinely tri-Hamiltonian system (it was known earlier to be a 
bi-Hamiltonian system). We have clarified the connection of this sytem to the 
modified dispersive water wave equation and have studied various properties 
of our model in the dispersionless limit. 
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